DoE 32: Analysing a 2^2 fatorial design in R